Atomistic study on the microscopic mechanism of grain boundary embrittlement 


induced by small dense helium bubbles in iron 


Lei Peng **, Yongjie Sun ^", Jingyi Shi^ ®*, Yifei Liu ^^, Shangming Chen ^, Liuliu Li ^^* 


^ State Key Laboratory of Particle Detection and Electronics, University of Science and 
Technology of China, Hefei, Anhui 230026, China 

^ School of Nuclear Science and Technology, University of Science and Technology of China, 
Hefei, Anhui 230027, China 

* Science and Technology on Reactor System Design Technology Laboratory, Nuclear Power 
Institute of China, Chengdu, Sichuan, 610213, China 


*Corresponding author. 
E-mail address: shijy@ustc.edu.cn (Jingyi Shi) 
This work was supported by the National Natural Science Foundation of China (Nos. 12175231 and 


11805131), Anhui Natural Science Foundation of China (No. 2108085J05), the National Key 
Research and Development Plan of China (No. 2018YFE0307101), and the Collaborative 
Innovation Program of the Hefei Science Center, CAS (No. 2021HSC-CIP020 and 2022HSC- 
CIP009). 


Abstract 

The helium bubbles induced by 14 MeV neutron irradiation can cause intergranular 
fractures in reduced activation ferritic martensitic (RAFM) steel, which is a candidate 
structural material for fusion reactors. In order to elucidate the susceptibility of different 
grain boundaries (GBs) to helium-induced embrittlement, the tensile fracture processes 
of 10 types of GBs with and without helium bubbles in body-centered cubic (bec) iron 
at the relevant service temperature of 600 K were investigated via molecular dynamics 
methods. The results indicate that in the absence of helium bubbles, the GBs studied 
here can be classified into two distinct categories: brittle GBs and ductile GBs. The 
atomic scale analysis shows that the plastic deformation of ductile GB at high 
temperatures originates from complex plastic deformation mechanisms, including the 
Bain/Burgers path phase transition and deformation twinning, in which the Bain path 
phase transition is the most dominant plastic deformation mechanism. However, the 
presence of helium bubbles severely inhibits the plastic deformation channels of the 
GBs, resulting in a significant decrease in elongation at fractures. For bubble-decorated 
GBs, the ultimate tensile strength increases with the increase of the misorientation angle. 
Interestingly, the coherent twin boundary X3(112j was found to maintain relatively 


high fracture strength and maximum failure strain under the influence of helium bubbles. 
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1. Introduction 


Reduced activation ferritic martensitic (RAFM) steel and its improved models, 
such as oxide dispersion strengthened (ODS), thermomechanical treatments (TMT), 
and castable nanostructured alloys (CNA) steels, are the primary candidate structural 
steels for future fusion reactor and advanced fission reactors owing to their excellent 
irradiation resistance, high temperature performance, and high level of technological 
and processing maturity [1—8]. However, in the service environment of fusion reactors, 
high-energy neutron irradiation will introduce a large number of helium atoms (i.e., up 
to 10-15 appm/dpa) into structural steel through the (n, a) reaction. [9-11]. Due to the 
extremely low solubility of helium atoms in metallic alloys, helium atoms have a strong 
tendency to coalesce in the form of nanometric bubbles. Moreover, helium bubbles tend 
to nucleate and grow at intrinsic or irradiation defects with sizeable free volumes, such 
as grain boundaries (GBs) and dislocations. While helium contributes to the safe use of 
fission reactor materials in some cases, such as by protecting the integrity of the fuel 
rod cladding [12], it is harmful to structural material used in fusion reactors. A high 
density of helium bubbles on the GB can significantly weaken the GB cohesive strength 
and lead to a drop in the fracture stress, exacerbating the failure of RAFM steel via 
intergranular fracture. Therefore, it is essential to elucidate the underlying mechanism 
of helium-induced GB embrittlement and establish the relationship between the helium 
embrittlement susceptibility of GB and GB structure characteristics to explore possible 
ways to improve the resistance of RAFM steel against helium-induced GB 
embrittlement. 

Few experiments have investigated the mechanical response of individual GBs to 
helium. For example, Miura et al. [13] investigated the effect of helium on the fracture 
properties of six different non-coincidence site lattice GBs of austenitic stainless steel 
by using micro-tensile tests combined with electron backscatter diffraction (EBSD) and 


focused ion beam (FIB) techniques. The results indicate that the segregation capacity 


2 


of GBs to helium impurities plays an important role in GB helium embrittlement. 
However, there is no report on the results of GB micro-tensile tests of RAFM steels. 
Atomic scale simulation methods are convenient for investigating the mechanical 
behavior of different types of GBs under the influence of helium bubbles. Terentyev et 
al. [14] used uniaxial tensile simulations to study the effect of helium bubbles on the 
fracture processes of six different GBs with «110» as the tilt axis at 0 K in body- 
centered cubic (bcc) iron. The results indicate that the helium bubbles have a substantial 
effect on the fracture stress and strain of GBs and an intensive suppression on the slip 
related to plastic deformation. Moreover, our previous work [15] also investigated the 
effect of small helium-vacancy clusters on the tensile deformation of two types of 
symmetric tilt grain boundaries (STGBs) with <110> and <100> tilt axes at both 0 K 
and 300 K. The results indicated that the helium-vacancy clusters exert a significant 
suppressive effect on GB sliding and that the size effect of helium-vacancy clusters on 
GB sliding seems to be stronger than that of the He/V ratio. However, whether the 
helium-embrittlement mechanism changes for different GBs remains unclear. 
Moreover, in these previous simulation studies, a single helium bubble was created 
by removing Fe atoms from a spherical region and filling up the volume with helium 
atoms, which did not take into account the effect of GBs on helium segregation. In our 
recent work, we comprehensively investigated the effects of GB characteristic 
parameters on helium segregation, helium bubble nucleation and growth processes in 
bcc iron for high-energy neutron irradiation-related service environments [16]. The 
spatial distributions of helium bubbles in bubble-decorated bicrystal models were 
obtained via kinetic evolution at 600 K, and the helium concentration was set as 2000 
appm. In this paper, the uniaxial tensile simulations of 10 different clean GBs (4 GBs 
with «100^ as the tilt axis and 6 GBs with «1107 as the tilt axis) and their corresponding 
bubble-decorated GBs were performed at 600 K using molecular dynamics methods. 
The experimental characterization of ferritic steels has shown that several of the studied 
GBs were observed at concentrations higher than those of random GBs [17, 18]. The 
microscopic mechanism of helium-induced GB embrittlement at high temperatures was 


studied by comparing and analyzing the tensile processes of clean GBs and bubble- 
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decorated GBs. The fracture strength and strain of bubble-decorated GBs were analyzed 
and correlated with the misorientation angle of GB to characterize the helium 


embrittlement susceptibilities of different GBs. 
2. Simulation methods 


In this work, the parallel molecular dynamics (MD) code LAMMPS [19, 20] was 
used to simulate the tensile processes of GBs in bcc iron with and without the presence 
of helium bubbles. The s-band Fe-He many-body potential [21], along with the Fe-Fe 
potential of Ackland and Mendelev [22] and He-He potential of Aziz [23, 24], was 
selected to describe the interatomic interaction in the Fe-He system. Based on the 
potential set, the calculated properties are in good agreement with data from 
experimental and ab initio calculations demonstrated in previous works [25—29]. 

Based on the coincidence site lattice (CSL) theory, 10 types of STGBs that were 
observed at concentrations higher than those of random GBs [17, 18] were constructed 
via the bicrystal method. The GB is located at the interface of two perfect grains with 
different orientations. The detailed parameters of these GBs are listed in Table 1. To 
acquire the equilibrium GB structures of the bicrystal system, we sampled over 10,000 
structures using rigid body translations. The displacement shift complete lattice was 
used to characterize the translational vector. The equilibrium structure was identified as 
the optimized lowest energy structure. Table 2 shows a comparison of the GB energy 
calculated in this paper with the results of DFT calculations [30—35] and other MD 
calculations [14, 36-38]. The majority of the GB energy values calculated in this work 
are consistent with other MD calculations but are slightly lower than the results 
obtained via DFT calculations. Details of the bicrystal method for the construction of 
GB structures can be found in Refs. [39-43]. In the Cartesian coordinate system of the 
simulation box, the tilt axis is aligned with the z-coordinate direction, and the normal 
direction of the GB is parallel to the y-coordinate direction. Fig. 1 shows the side view 
of the supercells for 10 different STGBs. The structural units that make up each GB are 


marked in the figure. 


Table 1. Detailed parameters of the studied STGBs. These parameters include tilt axis 
(<110> and <100>), coincidence index (2), GB plane index, misorientation angle (0) in 
degrees, total number of iron atoms (Nre), initial supercell dimensions (Lx, Ly, and Lz) 
in nanometer, number density of GB helium bubbles (Dune) in m? and average diameter 


of GB helium bubbles (D;) in nm. 


Tilt GB 0 Lx Ly Lz Due Da 
x Nee 
axis plane | (deg) (nm) | (nm) | (nm) | (1025 m?) | (107! nm) 
{111} | 109.47 | 300000 | 17.5 | 24.8 | 8.1 1.8 2.98 
3 
{112} | 70.53 | 302400 | 17.3 | 25.2 | 8.1 1.9 4.07 
{114} | 38.94 | 286400 | 17.1 | 24.2 | 8.1 2.7 3.53 
<110> | 9 
{221} | 141.06 | 300160 | 17.0 | 25.5 | 8.1 4.5 3.06 
{113} | 50.48 | 297440 | 17.4 | 24.7 | 8.1 2.7 3.57 
11 
{332} | 129.52 | 285120 | 17.0 | 24.1 | 8.1 2.5 3.94 
{210} | 53.13 | 299488 | 17.9 | 24.5 | 8.0 2.4 3.58 
5 
{310} | 36.87 | 297920 | 17.2 | 25.3 | 8.0 3.7 3.59 
<100> 
{320} | 67.38 | 296072 | 17.5 | 24.7 | 8.0 2.4 3.82 
13 
{510} | 22.62 | 313152 | 17.5 | 26.1 | 8.0 2.6 3.68 


Table 2. Comparison between the GB energies (Ecg, in J/m?) calculated in this work 


and in the literature. 


Ecs (J/m?) 


STGB type MD calculations DFT calculations 
X3«110» (111) 1.308 1.308", 1.295° 1.51°, 1.525, 1.538 
©3<110>{112} 0.260 0.260? , 0.262? 0.34*, 0.432, 0.45” 
x5«100» (210) 1.113 1.113? 1.618, 1.64" 
x5«100» (310) 1.008 1.008*, 1.19* 1.48*, 1.49f, 1.532 
£9<110>{114} 1.286 1.286%, 1.404 1.50%, 1.388 
£9<110>{221} 1.168 1.172°, 1.167” 1.66", 1.714, 1.62) 
£11<110>{113} 1.113 1.113, 1.03° 1.455, 1.44 
£11<110>{332} 1.019 1.022, 1.004 1.382, 1.475, 1.49! 
213«100» (320! 1.108 1.108? 1.53? 
£13<100>{510} 1.005 1.005? 1.40^ 


‘Ref. [36]; "Ref. [14]; *Ref. [37]; ‘Ref. [38]; “Ref. [30]; ‘Ref. [31]; Ref. [32]; "Ref. [33]; 
‘Ref, [34]; Ref. [35]. 
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Fig. 1. Schematic diagram of bcc iron bicrystals with STGBs located at the center. Iron atoms of 
each «100» or «110» layer are shown in a different color. The A, B, C and D notation refers to the 
different structural units that make up the GBs. The marks consisting of A, B, C and D below each 


GB refer to the structural units and its order that make up this GB. 


The He concentration in structural steels for the first wall of a fusion power plant 


can reach 1090 appm after 5 years [44] and 2000 appm after longer service times [10]. 
Thus, the initial helium concentration for our simulation was set as 2000 appm to study 
the GB mechanical properties at a high He concentration. Initially, all of the helium 
atoms were randomly distributed in the bicrystal model. A stable helium bubble 
distribution in the bicrystal model was obtained after 1.6 ns of MD evolution with the 
NPT ensemble under 600 K. Then, the bicrystal model with helium bubble distribution 
was used as the initial model for GB tensile simulation. The detailed process for 
building the initial models of bubble-decorated GBs and quantitative information on 
helium bubbles can be found in our previous work [16]. Owing to the high yield of 
helium, spallation neutron irradiation applications are widely used to test the properties 
of fusion reactor candidate materials. The distribution of high-density and small-sized 
helium bubbles on the GBs in the bubble-decorated GBs models (as shown in Fig. 2(b)) 
is consistent with the distribution of bubbles on the GBs in the spallation neutron 
irradiated RAFM steel specimens (as shown in Fig. 2(c, d)) presented in Refs. [45, 46]. 

The uniaxial tensile simulations for GBs with and without helium bubbles were 


! after 


performed under the same conditions. The tensile strain rate was set as 10° s 
comparison tests with 107 s! and 10? s'!. The results showed that the deformation 
mechanism of GBs remains unchanged at different strain rates. Under the strain rate of 
108 s^, relatively accurate fracture strength and strain values can be obtained while 
greatly saving computational resources. The direction of the uniaxial tensile load is 
perpendicular to the GB plane. Periodic boundary conditions were applied along x and 
z directions. To avoid the interactions of boundary atoms in the y-direction, iron atoms 
more than 120 nm away from the GB plane were removed to create a vacuum layer, as 
shown in Fig. 2. In order to control the tensile strain, iron atoms located 115—120 nm 
from the GB plane were fixed. The strain was increased by moving both of the fixed 
layers, one along the positive y-direction and the other along the negative y-direction, 
and remapping the rest of the atom coordinates. Each strain step increased the strain by 
0.01 percent, and this was followed by a 1.0 ps relaxation with the NVT ensemble. This 
process was executed cyclically until the bicrystal model fractured. After the relaxation 


of each strain step, the strain and stress in the y-direction were recorded to analyze the 
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tensile deformation process of the GB. The visualization and analysis of the atomic 
configurations were performed using the dislocation extraction analysis (DXA) and 


common neighbor analysis (CNA) as implemented in OVITO [47]. 
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Fig. 2. The uniaxial tensile model of (a) clean GB and (b) bubble-decorated GB, taking the 
%3<110>{112} GB as an example. The red spheres in panel (b) represent helium atoms, and the 
blue spheres in both figures represent iron atoms. The vacuum layers, fixed layers, and GB planes 
are indicated by black, red, and green rectangular boxes, respectively. Panels (c) and (d) [46] show 
a high-density distribution of small-scale helium bubbles at GBs in RAFM steel irradiated by 


spallation neutrons. 


3. Results 


The studied GBs can be classified into two categories based on the tensile behavior 
of clean GBs: (1) brittle GBs with no plastic deformation stage, including 231111], 
x9 {114}, and £9 {221}GB (hereafter, the notation of the tilt axis is omitted) and (2) 
ductile GBs with a plastic deformation stage. Ductile GBs have complex plastic 
deformation mechanisms at high temperatures. The plastic deformation mechanisms 
observed in the current work can be classified into three categories: (1) Bain path phase 
transition with a lath structure, including £5 {210}, &5 {310}, £13 {320} and 213 {510} 
GB; (2) deformation twinning, which only occurs during the uniaxial tensile process of 
the twin boundary %3{112}; (3) holonomic Bain/Burgers path phase transition, 
including X11 {113} and £11 {332} GB. In order to comprehensively study the plastic 
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deformation behaviors of different clean GBs and elucidate the influence of helium 
bubbles on tensile deformation behavior, the tensile processes of GBs were investigated 
in detail by combining the stress-strain curves with an atomic configuration analysis. 
The tensile deformation processes of X5 {210}, X5(310], 213 {320}, and X13 {510} 
GBs via the first deformation mechanism, i.e., Bain path phase transition with a lath 
structure, were analyzed in Fig. 3. The results indicate that during uniaxial tensile 
simulations, the four clean GBs all underwent a plastic deformation stage in which the 
stress first dropped rapidly and then increased slowly. £5 {310} GB was taken as an 
example to illustrate the deformation process from the atomic scale configuration 
perspective. Fig. 3(b-e) show the atomic scale configuration of the bicrystal at different 
stages during the tensile process. The stress and strain corresponding to each 
configuration are indicated with stars in Fig. 3(a). Fig. 3(b) shows that the initial GB 
configuration remains intact at high temperatures. In the elastic stage, the stress 
increased linearly as the strain increased, and no phase transition initially occurred. Fig. 
3(c) shows that once the strain reached 7.28% and the corresponding stress reached 
9.77 GPa, the GB region of the upper grain began to undergo a Bain path phase 
transition (bec to fcc). Fig. 3(d) shows that as the tensile strain increased, the phase 
transition region in the grain then began to expand to the fixed end in both upper and 
lower grains, and fcc lath grains were formed in the phase transition region. As the 
strain further increased, the phase transition region continued to expand, leading to an 
increased size of the lath grains. Fig. 3(1) shows the largest phase transition region 
covering almost the entire grains. Soon after the strain reached a limit of 18.95%, cracks 
started to form at the GB plane. However, Fig. 3(a) shows that with the presence of 
helium bubbles, there was no plastic deformation stage in any of the four GBs. The 
elongation at the fractures decreased sharply, and the ultimate tensile strength also 
decreased. Fig. 3(f-1) shows the atomic configurations of crack initiation in bubble- 
decorated £5 {210}, 251310], 213 {320}, and X131510) GB. As shown in the figure, 
no phase transition region appeared in the grains before the fracture of the bicrystals. 
When the applied stress surpassed a certain threshold, the dislocation loops nucleated 


around the helium bubbles in the bulk gradually grew and finally reached the fixed end 
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of the bicrystal model, which absorbed part of the strain energy. Moreover, the helium 
bubbles on the GB reduced the bonding area of the boundary, which led to a drop in 
cohesive strength. In this case, the fracture stress was lower than the critical stress for 
the Bain path phase transition. This means that the existence of helium bubbles at the 
GB inhibited the channel of the Bain path phase transition, which led to the brittle 


rupture of the bicrystal at the GB. 
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Fig. 3. Stress-strain curves and atomic configurations of GBs with a Bain path phase transition and 
lath structure. (a) Stress-strain curves of £5 {210}, £5 {310}, 213 {320}, and £13 {510} GBs with 
(dotted line) and without (solid line) helium bubbles. Panels (b-e) illustrate the atomic 
configurations of clean £5 {310} GB at different deformation stages. The blue, green, red, and white 
spheres represent iron atoms with bcc, fcc, hcp, and other structures, respectively. Panels (f-i) 
illustrate the atomic configurations of the crack initiation of bubble-decorated 25 {210}, £5 {310}, 
X13 {320}, and X13 {510} GBs, respectively. Parts of the curves are horizontally translated for clarity. 
The blue spheres, red spheres, green lines, and pink lines represent iron atoms, helium atoms, 2<111> 
dislocation/dislocation loops, and <100> dislocation/dislocation loops, respectively. The black 


rectangle indicates the regions where crack initiation occurred. 


The tensile deformation process of X3[112) GB via the second deformation 
mechanism, i.e., deformation twinning, was analyzed, as shown in Fig. 4. X3 {112} GB 
is a unique twin GB, and the tensile direction along the <112> orientation leads to the 
activation of its twin system {112}<111> [48, 49]. In particular, Fig. 4(a) depicts the 
stress-strain curve of clean £3 {112} GB, which shows that the tensile stress presented 


a fluctuating increase at the plastic deformation stage. During the uniaxial tensile test, 
10 


atomic configurations of clean £3 {112} GB were analyzed to explain this phenomenon, 
as shown in Fig. 4(b—f). The stress and strain corresponding to each configuration are 
indicated in Fig. 4(a) by stars. In contrast to other GBs, the initial £3{112} GB 
configuration is a narrow region consisting of two atomic layers, which reflects the 
peculiarity of &3{112} GB as a twin GB, as shown in Fig. 4(b). The stress increased 
linearly with the strain in the elastic stage, and no phase transition occurred at this stage. 
Fig. 4(c, d) show that when the tensile strain reached 8.13% and the corresponding 
stress rose to 12.78 GPa, the GB plane started to slip along the (1-21)[111] slip system 
to the fixed end and formed a set of deformation twins perpendicular to the original GB 
plane. The GB structures of the deformation twins were entirely consistent with the 
original twin boundary. Thereafter, the shape and size of the deformation twins were 
gradually optimized to be consistent under the action of applied stress, as shown in Fig. 
4(e). Fig. 4(h) shows that the structures of these deformation twins closely match with 
those observed in RAFM [50, 51] and austenitic steels [13], as observed by TEM. 
However, the difference in critical stress for deformation twinning was relatively large, 
and this was mainly owing to the high strain rate in the MD simulation. During the 
twinning period, the tensile stress rose in a fluctuating manner. Fig. 4(f) shows that with 
the further increase of strain, the high-stress state at the fixed end led to the formation 
of a Bain path phase transition region near the newly formed twin boundaries. 
Deformation twinning refined the newly formed grains, quickly hindering the 
expansion of the phase transition region. Eventually, the bicrystal fractures at the 
interface between the deformation twin boundary and the fixed end, which is consistent 
with the experimental finding that deformation twins can act as micro-crack initiators, 
as shown in Fig. 4(1) [51]. The stress-strain curve of the bubble-decorated £3 {112} GB 
is shown with a dashed line in Fig. 4(a). However, compared with clean £3 {112} GB, 
the fracture strength and elongation of bubble-decorated GB are much lower. The 
fracture stress is lower than the critical stress for deformation twinning. Fig. 4(a) shows 
that despite the process of deformation twinning having been inhibited, the twin 
boundary still retained some plasticity under the influence of helium bubbles. The 


atomic-scale configuration showed that the plasticity derived from the relatively active 
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dislocation evolution process at the crack tip, as shown in Fig. 4(g). 
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Fig. 4. Stress-strain curves and atomic configurations of 331112] GB with deformation twins. (a) 
Stress-strain curves of 331112) GBs with (dotted line) and without (solid line) helium bubbles. 
Panels (b—f) illustrate the atomic configurations of clean £3 {112} GB at different deformation stages. 
The blue, green, red, and white spheres represent iron atoms with bcc, fcc, hcp, and other structures, 
respectively. Panel (g) illustrates the atomic configuration of the crack initiation of helium bubble- 
decorated £3{112} GB. The blue spheres, red spheres, green lines, and pink lines represent iron 
atoms, helium atoms, %2<111> dislocation/dislocation loops, and <100> dislocation/dislocation 
loops, respectively. Panel (h) [51] shows the deformation twin structure observed during tensile tests 
under TEM. Panel (1) [51] shows the initiation of microcracks caused by deformation twins. 


Fig. 5 shows the tensile deformation processes of £11{113} and £11 {332} GBs 
via the third deformation mechanism, i.e., holonomic Bain/Burgers path phase 
transition. Fig. 5(a) shows that both clean GBs underwent a plastic deformation stage 
with constant stress during the uniaxial tensile simulations. However, the atomic scale 
analysis revealed that the two GBs had different phase transition mechanisms. Fig. 5(b— 
e) show the atomic scale configuration of the X11 {332} GB with a Burgers path phase 
transition (bcc to hcp) at different stages during the tensile process. In Fig. 5(a), a star 
indicates the corresponding stress and strain of each configuration. Fig. 5(g—j) show the 
atomic scale configuration of the £11{113} GB with a Bain path phase transition at 


different stages during the tensile process. The corresponding stress and strain of each 
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configuration are also indicated in Fig. 5(a) by stars. As shown in Fig. 5(b, g), the initial 
configurations of X11(332) and X11{113} GB were different. Some atoms at the 
X11 {332} GB plane had atomic configurations similar to those of hcp crystals, which 
provided a large number of nucleation sites for Burgers path phase transitions. 
Nonetheless, the £11 {113} GB does not have such a structure. As shown in Fig. 5(c, h), 
when the strain increased to 7.64%/7.0%, the Burgers/Bain path phase transition began 
to progress steadily in the direction perpendicular to the original GB plane in X11 {332} 
and X111113) GB, respectively. As strain increased, atoms with the bcc structure on 
both sides of the GB continuously transformed into the hcp or fcc structure. 
Simultaneously, the bcc-hcp or bcc-fcc phase boundary began to move toward the fixed 
end of the bicrystal. During this period, the stress remained almost constant as the strain 
increased. Then, once the strain reached 10.86%/13.00%, the volume of the hcp/fcc 
region reached its maximum value for X111332] and X111(113) GB, respectively, as 
shown in Fig. 5(d, 1). 

Typically, the phase transition process continues until the bcc-hcp or bcc-fcc phase 
boundary reaches the fixed end at 0 K [15]. However, the phase boundary is unstable 
at high temperatures. Fig. 5(e) shows that for X11(332]) GB, a stacking fault slip 
occurred at the upper bec-hcp phase boundary, resulting in the disruption of the Burgers 
path phase transition and the further destruction of the phase transition region. Likewise, 
Fig. 5(j) shows that for £11{113} GB, a series of dislocations were emitted from the 
lower bcc-fcc phase boundary and piled up at the fixed end of the lower grain, further 
resulting in the stress concentration and grain cracking at the lower fixed end. Also, 
under the effect of high temperature, at the beginning of the plastic deformation of the 
X111113j GB, there was a stage where the stress briefly decreased. This stage was 
related to one dislocation emission at the GB in the initial stage. After the dislocation 
emission, the phase transition region was divided into two parts, as shown in Fig. 5(1). 
Finally, when the strain reached 10.96%/13.42%, an initial crack formed in the slip band 
for X11 {332} GB or at the lower fixed end for X111(113) GB, as shown in Fig. 5(e, j). 

In contrast, there was no phase transition during the tensile process of bubble- 


decorated X11(332j and X111(113j GB, and the initial crack presented at the grain 
13 


boundary, as shown in Fig. 5(f, k). Likewise, the fracture stress was lower than the 
critical stress for the Bain/Burgers path phase transitions, which led to the brittle 


fracture and a noticeable drop in the elongation. 
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Fig. 5. Stress-strain curves and atomic configurations of GBs with holonomic Bain/Burgers path 
phase transitions. (a) Stress-strain curves of £11 {113} and X111(332) GBs with (dotted line) and 
without (solid line) helium bubbles. Panels (b-e) and (g-j) illustrate the atomic configurations of 
clean X11(332] and £11{113} GBs at different deformation stages, respectively. The blue, red, 
green, and white spheres represent iron atoms with bcc, hcp, fcc, and other structures, respectively. 
Panels (f) and (k) illustrate the atomic configurations of the crack initiation of helium bubble- 
decorated GBs. The blue spheres, red spheres, and green lines represent iron atoms, helium atoms, 
and 1/2<111> dislocation/dislocation loops, respectively. 


Fig. 6 shows the fracture processes of the three brittle GBs, namely, X31(111j, 
x9 {114}, and X91221], with and without helium bubbles. Fig. 6(a) shows that the three 
clean GBs have similar stress-strain curves. Fig. 6(b-d) present the atomic 
configurations of three clean GBs just before fracture. The stress and strain 
corresponding to each configuration are indicated in Fig. 6(a) by stars. The 
configurations show that under the combined effects of high temperatures and stress, 
the GB regions of these bicrystals widened, and their original GB structure units were 
completely lost. In fact, the spatial distribution of atoms in the GB region became 
abnormally disordered [16]. The disorder of the GB structure hindered the 
aforementioned phase transition paths. Therefore, for clean GBs, the disorder of GB 
atoms is the fundamental cause of their brittle fractures. Fig. 6(e, f) show that the tensile 


behavior of bubble-decorated brittle GBs is similar to that of bubble-decorated ductile 
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GB. Under the influence of helium bubbles, the fracture strength and elongation of GBs 


both decreased to varying degrees, as expected. 
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Fig. 6. Stress-strain curves and atomic configurations of brittle GBs. (a) Stress-strain curves of 
x3 {111}, 291114], and £9{221} GBs with (dotted line) and without (solid line) helium bubbles. 
Panels (b-d) illustrate the atomic configurations of clean £3{111}, 29{114}, and X9(221) GBs 
under ultimate tensile stress, respectively. The blue, green, and white spheres represent iron atoms 
with bcc, fcc, and other structures, respectively. Panels (e—g) present the atomic configurations of 
the crack initiation of bubble-decorated 231111], X9(114), and £9{221} GBs, respectively. The 
blue spheres, red spheres, and green lines represent iron atoms, helium atoms, and 1/2<111> 


dislocation/dislocation loops, respectively. 


4. Discussion 


To further analyze the role of the Bain path phase transition in the plastic 
deformation of ductile GBs, the orientation relationship in the Bain path phase 
transition from bcc to fcc was investigated, as shown in Fig. 7. Fig. 7(a—e) shows the 
correspondence of the crystal planes and crystallographic directions in the bcc and fcc 
phases for 35(210], X5{310}, X13(320], X131510], and £11{113} GBs during the 
phase transition, respectively. In the phase transition process from the bcc structure to 
fcc structure, the «110», direction in the bcc structure and the «121-7, direction in the 
fcc structure were always parallel. In addition, the {110}. plane in the bcc structure was 


always parallel to the {111}, plane in the fcc structure. Fig. 7(f) presents a schematic 
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diagram of the Bain path phase transition and shows that as long as a tensile 
deformation of about 40% occurs along the y-axis (<100> crystal direction), that is, as 
long as the lattice constant in the y-direction is elongated from the original a to V2a, 
the original bcc unit cell can be transformed into an fcc unit cell. The newly formed fcc 
unit cell has a lattice constant of 2a, and can therefore store higher strain energy than 
the original bcc unit cell with a as its lattice constant. For crystal planes, the {111}, 
planes in the fcc lattice corresponded to the {110}, planes in the bcc lattice. For the 
crystal orientation, the «112», orientation in the fcc lattice corresponded to the «110»; 
orientation in the bcc lattice, and the <110>, orientation in the fcc lattice corresponded 
to the <111>, orientation in the bcc lattice. The correspondence of the crystal planes 
and crystallographic directions in the bcc and fcc phases during the GB phase transition 


are in agreement with the Bain path phase transition mechanism found in the bulk [52]. 
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Fig. 7. Schematic diagram of the orientation relationship in the Bain path phase transition from a 
bcc to an fcc lattice. (a) X5(210); (b) X5(310); (c) £13 {320}; (d) X13(5103; (e) X1111133; (f) 
correspondence of the crystal planes and crystallographic directions between the bcc and fcc phases. 
The blue, green, and yellow spheres represent iron atoms with bcc, fcc, and hcp structures, 
respectively, in panels (a-e). The solid red spheres represent the original bcc lattice before the phase 


transition in panel (f). 


In order to investigate the susceptibility of different GBs to helium embrittlement, 
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the relationship between the ultimate tensile (fracture) strength as well as the failure 
(fracture) strain and the misorientation angle under the influence of helium bubbles was 
analyzed, as shown in Fig. 8. For a comparative analysis, the ultimate tensile strength 
and failure strain for clean GBs are also included in Fig. 8(a, b). First, the failure strain 
dropped under the influence of helium bubbles. Second, the ultimate tensile strength of 
individual GBs also decreased under the influence of helium bubbles, and this finding 
has been confirmed by the results of Ref. [13]. Previous experimentation has shown 
that the ultimate tensile strength of helium-containing GBs is related to the 
misorientation angle [13]. Our results in Fig. 8(a) further show that the ultimate tensile 
strength generally displayed an increasing trend as the misorientation angle increased. 
There is an excellent linear relationship between the ultimate tensile strength and the 
misorientation angle, especially for GBs with a «100» tilt axis. Compared with the 
tensile results of clean GBs, it was found that although helium bubbles reduce the 
ultimate tensile strength of each GB, they do not break the relationship between the 


ultimate tensile strength and misorientation angle. 


Fig. 8(b) shows the relationship between the failure strain and misorientation angle. 
The coherent twin boundary 231112], with its low GB energy and sigma value, had a 
larger failure strain than did the other bubble-decorated GBs, and it also had the lowest 
strain drop among the ductile GBs. Apart from the £3 {112} GB, the fracture strain 
generally decreased as the misorientation angle increased. Compared with other bubble- 
decorated GBs, 231112) GB displayed a high ultimate tensile strength while 
maintaining the largest failure strain. Therefore, it is suggested that the distribution of 
x3 {112} GBs in RAFM steel should be optimized by grain boundary engineering to 
improve their resistance to GB helium embrittlement, similar to the related 


improvement of resistance against intergranular damage in Refs. [8, 53—55]. 
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Fig. 8. The relationship between (a) the ultimate tensile strength and (b) the failure strain and the 
misorientation angle for clean (hollow) and bubble-decorated (solid) GBs. The red circles and green 
squares refer to <110> and <100> tilt axis GBs, respectively. 


5. Conclusions 


In order to elucidate the microscopic mechanism of high temperature GB helium 
embrittlement and the susceptibility of different types of GBs to helium embrittlement 
in RAFM steel, we investigated the uniaxial tensile processes of 10 different clean GBs 
with two types of tilt axis (1.e., «1107 and <100>) as well as their corresponding bubble- 
decorated GBs in bcc iron at 600 K. The tensile behavior reveals that the studied GBs 
fall into two distinct categories: brittle and ductile. There are three different types of 
plastic deformation mechanisms in ductile GBs. The first type is a Bain path phase 
transition with a lath structure, which only occurs during the tensile process of GBs 
with <100> as the tilt axis. The second type is deformation twinning, which only occurs 
at the unique twin GB &3 {112}. The third type is a holonomic Bain/Burgers path phase 
transition, which occurs in some of the GBs with <110> as the tilt axis. The results ofa 
structural analysis of the tensile process at the atomic scale imply that the plastic 
deformation mechanisms of almost all of the ductile GBs studied here are related to the 
Bain path phase transition, which plays a vital role in the plastic deformation of the 
GBs. However, the presence of helium bubbles results in all GBs failing via brittle 
fracture and causes both the GB ultimate tensile strength and the failure strain to 


decrease. Helium bubbles significantly reduce the failure strain of ductile GBs by 
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inhibiting the channel of plastic deformation. In the presence of helium bubbles, the 
ultimate tensile strength of GBs increases with the increase of the GB misorientation 
angle. The coherent twin GB 231112) maintains the maximum failure strain and a 


relatively high ultimate tensile strength under the influence of helium bubbles. 
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